AD-A097  616  ROT/U.  AIRCRAFT  ESTABLISHMENT  FARNB0R0U6H  (EN6LANOI  F/6  1/1 

NOVL60ICVLF,SoLCHa5oVION  **  **  C«**»”ERISTICS  OF  A  -CTC 

UNCLASSIFIED  RAC -LIBRARY  TRANS-2060  DRIC-M-7T4M 


DTIC 


Trans  2060 


UNLIMITED 


tfe 

W  #  r 

ROYAL  AIRCRAFT  ESTABLISHMENT 


B ii776': 0  - 

Trans  2060 


a> 


Library  Translation  2060 

November  1980 

ANALYTICAL  INVESTIGATION  OF  THE 
NONLINEAR  CHARACTERISTICS  OF  A 
RECTANGULAR  WING  OF 
SMALL  ASPECT  RATIO 


V.F.  Molchanov 


LCf  VviG.IT  t;  1 

PJ!  I 

i 

rr\.- ' ; 

^  •  -  v  .  :a  I 


<8  i  3 


Procurement  Executive,  Ministry  of  Defence 
Farnborough,  Hants 


rN  -  , 

f 


i  v>  r 


UNLIMITED 


Translations  in  this  series  are  available 
from: 

THE  R.A.E.  LIBRARY 
Q.4  BUILDING 

R.A.E.  FARNBOROUGH 
HANTS 

New  translations  are  announced  monthly  in 


"LIST  OF  R.A.E.  TECHNICAL  REPORTS, 
TRANSLATIONS  and  BIBLIOGRAPHIES" 


\ 


A k  'LI.  i  r  A  h  /  T /,  ////.  -  .' .// / 


/  Hi 


1‘  1)1  S  M  .  h‘H  .  h 


K_;t_Y  \  1. _ ,\  I  R  C  R  A  ,•  [  H  S  T  \  B  I.  1  SUM  L  N  T 

l.ibrarv  Translation  2060 

.  Ilj 

Received  lor  nrintiru*  Id  November  -1981/ 


(:) 


analytical  investigation  of  the  nonlinear  characteristics  of  a 

JICOIANGULAR  TING  OF  SMALL  ASPECT  RATIO 

l  ANAL  I  T I CHKSKOK  „I  SSLEUOVAN I  E  ME  1.1  NIC  1 NYKH  KHARAK.TKRLST  I  K 
PRYAMOUCOL’NOGO  KRYLA  HALOGO  COL  I.VEX I YA)  » 


\b 


A  .a  ,  -r  V.  F. /loir 

<>  1  Ifton ii  ' 


by 

■lolrhanov 


1  r  JT’ii 1  tW 


Sarbarn  Cross  land  # 


,V,  IX,  No.  S,  1-10  (  1978) 

cj 

px-l.y 


Translation  editor 
.  1 .  1 1 .  H .  Smi  t  li 


AUTHOR' S  SUMMARY 

In  Ref'  I,  the  form  was  found  of  the  main  nonlinear  term  of  the  expansion  of  the 
lift  eoef f i eient  for  a  rectangular  wing.  In  this  paper,  the  form  is  found  of  all  terms 
>f  the  expansion  of  the  lift  and  moment  coefficients.  Generalisations  are  i>iven  for  the 
ase  of  certain  non-steady  flows.  Results  of  calculations  are  presented. 
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I  Non  1  i  ilea  r  effects  of  forces  ot  in  ideal  medium  .  n  t  i  n  >•.  on  ,1  wing  am,  in  genera  1 , 

governed  by  the  presenn'  at  a  surface  at  >1  i  scon  t  i  tin  i  t  y  at  t  tin  t  .infant  i  a  I  components  of 

> 

tlia  velocity.  Within  an  accuracy  of  the  order  at  >•  ,  these  effects  can  be  found  by  the 

method  of  plane  sections.  Here  t  is  the  angle  ol  incidence,  is  the  aspect  ratio, 

i  varies  within  the  limits: 


where  is  a  constant. 

i  ■> 

By  virtue  of  (l-l),  o'"  ^  t"  ,  which  makes  it  possible,  in  calculation,  to  replace 

sin  i  by  a  .  The  problem  cont'erniit,;  the  movement  of  a  rectangular  wing  comes  down  to 
one  of  finding  a  plane,  non-steady  separated  flow  of  an  ideal  liquid  about  a  plate 
instantaneously  brought  into  motion  along  its  normal,  with  a  velocity  v^  .  Thereupon, 
the  following  relationship  arises: 


iV  =  v  ,  tV  =  It 

n 


(1-2) 


where  v  is  the  velocity  of  motion  of  the  plate;  h  is  the  distance  between  the  leading 
edge  of  the  wing  and  the  plane  of  the  section  of  the  wing  normal  to  the  flow;  V  is  the 
velocity  at  infinity;  t  is  the  time  elapsed  from  the  moment  of  the  start  of  the  motion 
of  particles  of  the  medium  in  the  given  section. 


The  case  of  non- steady  movement  of  a  rectangular  wing  also  reduces  to  this  plane 
problem  when  the  wing  is  momentarily  brought  into  motion  with  a  constant  velocity  V  . 

If  t  fl  is  the  time  elapsed  from  the  start  of  the  wing  motion,  and  is  the  chord 

of  t /,e  wing,  then: 


for  0  <  h  <  V  t 


0 


(1-3) 


f  or  V  _  t  <  h  <  h 


0 


Then,  as  before,  t  is  the  time  elapsed  from  the  start  of  movement  of  particles  of  the 
medium  in  the  given  section. 


If  the  wing  is  placed  in  a  shock  tube,  and  if  t  is  the  lime  elapsed  from  the 
start  of  movement  of  the  particles  of  the  medium  near  the  leading  edge  of  the  wing,  then 
for  the  sections: 


0  <  h  <  V  t„ 

TO  0 


h 

V 


(  1-4) 


while,  in  sections  V_t^  h  <  h^  ,  motion  is  absent.  In  fact,  in  the  vicinity  of 
h  =  V  t^  ,  there  is  a  transition  zone.  However,  its  rejection  leads  to  an  error  of  the 
same  order,  . 

The  solution  of  the  two-dimensional  problem  makes  it  possible  to  find  the  drag  force 
on  the  plate  Y(t )  .  This  force  is  numerical lv  equal  to  the  normal  force  acting  on  a 


1 


small  part  of  t  lit'  wing  included  between  sections  h  and  h  +  h  ,  relative  to  the 
distance  between  the  sections  .'li  .  Therefore,  the  following  expressions  are  obtained  lor 
the  lift  coefficient  c  and  for  the  moment  coefficient  M  ,  relative  to  the  leading 

y  7. 

edge . 


0 


72  i 


Ydh  , 


■  v  S 


0 

I 


liYdh  , 


(  1-5) 


where  Y  must  be  expressed  in  terms  of  h  ,  V  and  t()  . 


In  practice,  it  is  more  convenient  to  use  the  impulse  I)  <  t )  instead  ol  the  force 


Y 1 1 ) 


D  (  ) 


Yd  t  . 


(  1-6) 


Thus  the  problem  of  investigation  of  the  characteristics  of  the  wing  c  and  reduces 

to  the  determination  of  the  impulse  Oft)  as  a  function  of  the  time,  for  the  two- 
dimensional  ease  of  non-steady  separated  flow  of  a  plane  flow  of  an  ideal  liquid  around  a 
plate. 

2  If  the  surfaces  of  tangential  discontinuity  formed  near  the  edge  of  the  plate  are 

known,  then  the  velocit;  field  can  be  found  according  to  known  formulae.  Therefore,  the 
problem  of  finding  D(t)  comes  down  to  finding  these  surfaces. 

The  equations  of  motion  for  surfaces  of  tangential  discontinuity  can  be  written  in 
the  following  l .  rm: 


L'C  ,t)  +  <vr\t)> 


dw* 

"dt 


(2-1) 


w<r,t)  ,  rQ  >  r  >  o  ,  rQ(t) 


imfufr  ,t)  +  <  v(r  ,o> )  =  o 


(2-2) 


Here  W  is  the  point  of  the  complex  plane  lying  at  the  surface  of  the  discontinuity;  T 
is  the  potential  jump  at  point  W  ;  T^  is  the  potential  jump  at  the  edge  of  the  plate; 

T  ’ ;  the  complex  velocity  of  the  non-separated  flow  of  the  stream  v  /i  about  the  plate, 
directed  towards  the  plate  along  the  normal;  V  is  the  complex  velocity  induced  by 
discontinuities  in  the  presence  of  the  plate  (discontinuities  adjoin  its  edges);  <V>  is 
the  half-sum  of  the  values  V  ,  calculated  for  both  sides  of  the  di sconl i nui t v  at  the 
point  W  .  The  superior  asterisk  denotes  a  complex-conjugate  value. 

The  derivative  in  the  right-hand  part  of  (2-1)  applies  for  the  condition  2  =  const. 

Equation  (2-1)  is  obtained  on  the  basis  of  the  known  properties  of  a  tangential  discontin¬ 

uity,  its  impermeability  and  the  absence  of  a  pressure  jump.  Equation  (2-2)  is  obtained 
on  the  basis  of  Zhukovskii's  postulate.  The  systems  (2-1)  and  (2-2)  are  complete  in  that 

they  contain  the  functions  W(T,t)  and  Tg(t)  •  but  they  must  be  supplemented  by  a 


I.T  .’060 


\ 
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relationship  between  W(’,t)  and  the  functions  T  and  (V)  .  This  relationship  is 
conveniently  found  by  introducing  tin*  auxiliary  complex  plane  z  .  Without  destroying  the 
generality,  it  is  possible  to  consider  the  half-width  of  the  wing  as  being  equal  to  unity, 
and  the  edges  as  lying  at  points  1+  and  -1  in  the  complex  plane  W  .  Then  z  is 
defined  bv  the  function: 

/  2 

W  =  .1  +  z  ,  (2-  i) 


mapping  the  exterior  of  the  segment  (-1,  +1)  of  the  plane  z  on  the  exterior  of  the 
segment  ('  —  i ,  +i)  of  the  plane  W  . 

The  relationship  between  W(T,t)  and  l’,  W(T,t)  and  (V)  is  given  by  the  follow¬ 
ing  expressions: 

l'  '  “S  *  <V>  =  <v>^  >  (2“4) 


u 


v 

n 

i 


0 


0 


1 

2-u 


(2-5) 


Here  the  inferior  asterisk  indicates  a  dependence  of  the  given  function  on  the  variable  of 
i n t egra t ion . 

The  expression  for  D(t)  is  obtained  on  the  basis  of  formulae  in  Ref  2: 


D  =  ttv  +  AD  , 
n 

r 

'o 

AD  =  2pRe  j  z*dr 
U 


(2-6) 


(2-7) 


Let  us  change  to  new  variables.  Let  us  suppose  that: 


0  <  w  <  1  .  (2-8) 


Moreover,  in  order  to  abbreviate  the  writing  of  functions,  !'  and  z  will  be  given  a 
supplementary  definition  for  the  ease  of  negative  values  of  •.  ,  in  accordance  with  the 
following  expressions: 


V(-u)  =  "(w)  , 


z(-w) 


* 

-Z  (  ,  ) 


(2-9) 


In  terms  of  these  variables,  the  equations  of  motion  take  the  following  form: 


V 


(p(0)  -  p  ( z. ) )  K  +  (1  -  .  )  — -  -  p ( 0)  =  0  , 


+  1 


p(z) 


p(0) 


-I 


J 

-1 


F.  = 


>  (2-10) 


i  +  ; r  i 

—  i 

z  : 


,  p  ( o ) 


(2-1  I ) 


i/  .  Re  I 

"  i 

-I 


z*cl'-  • 


(2-12) 


where  ft.  =  .‘(w)  =  1  for  w  >  0  and  -I  for  w  <  0  .  The  problem  reduces  to  the 
solution  of  system  (2-10)  and  the  introduction  of  results  in  (2-11)  and  (2-12). 

3  Let  us  investigate  the  principle  of  finding  the  solution. 

Let  the  solution  of  a  certain  physical  problem  reduce  to  the  solution  of  the 
operator  equation: 

Fz  =  0  ,  z  =  z(w, s)  ,  ( 3- 1  ) 


which,  by  assumption,  has  a  unique  solution,  and  this  solution  can  be  represented  in  the 
following  way: 


z 


Ea  :  ,  a  (w),  f  (s) 

n  n  n  n 

n=0 


(3-2) 


(.’ ' :  'rr'.i  kc  t  -.'nof  -'ow:c:  icr  haa  /•<„•.  •>;  m.k  >: 


$  0 
n 


( 3-3) 


;  *  0  for  s  —  0 

n 


(  3-4) 


n+  I 


f  or 


s  -  C 


n  =  0,1,2,...  .  (3-3) 


It  is  necessary  to  find  the  functions  $  (s)  .  In  the  general  case,  for  this  purpose 
equation  (3-1)  must  be  solved.  However,  if  the  operator  F  has  certain  special  proper¬ 
ties  and  if,  from  physical  considerations,  one  can  obtain  additional  restrictions  on  the 

desired  solution  z  ,  then  the  function  i)1  can  be  found  even  without  the  direct  solution 

n 

of  (3-1 ) . 

Let  us  consider  the  case  when  the  operator  F  has  the  following  three  properties. 


(I)  For  any  function  of  z  which  can  be  represented  as  in  (3-2)  to  (3-5)  and  which 
belongs  to  a  certain  acceptable  region  of  definition,  the  expression  Fz  can  be  written 
in  the  following  form: 


IT  3060 
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Fz  = 


)  (A.i)C  : 

t_,  11  11 

n=0 


<  i-fi> 


here  a  and  ;  are  infinitely-dimensioned  vectors  with  the  coordinates  a  and  ; 

'  n  n 

respectively;  A  and  :■  are  certain  operators  acting  on  a  and  ;  .  We  observe  that, 
using  only  the  condition  (1-5),  it  is  not  possible  to  place  the  expressions  I  •(  as  a 
function  of  s  in  increasing  order.  It  is  noL,  for  example,  possible  to  compare  tilt- 
orders  of  the  functions  and  ;  d$  /ds  .  However,  in  expression  ('1-6)  there  are  also 

terms  which  can  be  arranged  in  an  order  relationship  established  on  the  basis  of 
conditions  (3-5). 

(2)  If,  using  conditions  (3-5),  any  term  of  (3-6)  having  a  higher  order  than  some  other 
term  in  the  latter  expression  is  rejected,  then  only  a  finite  number  of  terms  depending 

remain . 


solely  on  a^  and 


'0 


For  the  sake  of  simplicity,  let  us  consider  that  only  two  terms  remain.  Let  us  call 
this  residue  the  dominant  part  (DP): 


DP(Fz)  =  (A|a0)(f|'t>0)  +  (A2a0)Cy!>0) 


(3-7) 


>:■:  t‘  %crr- 


(3)  The  dominant  part  of  the  increment  F z  -  F(z  -  zR)  ,  where: 


f  ..  , 

Lj  n  n 

n=k 


(3-8) 


also  consists  of  a  finite  number  of  terms,  each  of  which  depends  only  on  a^  ,  aR  ,  , 

:  .  Here  the  dependence  on  a  ,  t  is  linear,  and  the  operators  determining  this 

K  K  K 

dependence  are  independent  of  k  .  Accordingly,  we  shall  investigate  the  case  when  the 
dominant  part  of  the  increment  consists  only  of  two  terms: 

DPfFz  -  F ( z  -  zR)j  =  (a|  (a0)ak^:  i  (  V°k)  +  (A2(a0)akH:,2(‘;0)<1,k^  4  (3_9) 

Depending  on  the  vector  ,  in  the  expression  (3-7),  the  dominant  term  will  be  either 
the  first,  or  the  second,  or  both  terms  will  have  the  same  orders.  The  expression  (3-6) 
may  similarly  equal  zero  only  in  the  case  when  the  sum  of  each  one  of  the  dominant  terms 
is  separately  equal  to  zero.  Therefore,  can  satisfy  one  of  the  following  equations: 


A|a0 


A2a0 


(3-10) 


0 


<  3-  I  I  ) 


0 


8 


('V'o^’i  +  A2;i( 


(W/(  ,'2'!0) 


I  ’ 


=  COIl.S  t  . 


(3-12) 


It  is  possible  to  give  preference  to  one  of  the  equations  (3-10),  (  i- 1  I )  or  ( 3—12) 
on  the  basis  of  aildition.il  restrictions  on  the  solution  z  obtained  from  physical  consid¬ 
erations.  In  the  given  case,  the  condition  of  a  unique  solution  is  sufficient  to  justify 
rejection  of  equations  (  l-M'l  and  (5-11),  since  they  do  not  determine  the  t  unction  :  . 

In  the  third  case,  however,  a  unique  solution  is  possible. 

[f  (5-7)  were  to  consist  of  n  terms,  then  it  would  he  necessary  to  investigate 
Jn  -  I  equations.  Functions  a  ,  f  are  found  in  a  similar  manner.  It  is  sufficient 
merely  to  carry  out  a  substitution  of  the  variables  z.  «  ,i  +  r.  ,  where  z  is 

determined  by  the  sum  0-8),  and  to  require  that  the  dominant  member  of  the  expression 
Fta^r^  +  z.j)  =  0  .  Thus: 


F  ( a 


O’O 


+  z ,  )  = 


Ha0*V 


IF  (a 


0V0 


V  -  F(aoV- 


0-13) 


and,  as  soon  as  a^  and  f  are  known,  then  the  dominant  member  (DM)  of  the  first  term 
can  be  picked  out  directly.  Moreover,  by  virtue  of  (3-12),  its  order  is  higher  than  zero. 
Taking  condition  (3-6)  into  consideration,  let  us  suppose  that: 


DMfF(.Vo))  =  . 

From  the  expressions  contained  in  brackets  in  (3—1 3 )  we  preferentially  take  out  the 
dominant  part.  Taking  (3-9)  into  consideration,  we  obtain: 


DP  f  F  (; 


0  0 


z,)  -  F(a()f0)j 


iA;<vvO;(vh)  +  li-ir,) 


If  aQ-g  is  no1  the  solution  of  problem  i  5-1),  then  Bj(a^)  ^  0  and  consequently,  the 
order  of  (3-14)  cannot  be  less  than  the  order  of  the  terms  in  (3-15).  On  the  other  hand, 
none  of  the  terms  in  (  5-15)  can  have  an  order  less  than  that  of  (3-14),  since,  in  such  a 
case,  for  a  the  trivial  solution  ,i  ^  7  0  would  he  found.  One  of  the  following  three 
systems  of  equations  is  possible  for  a  : 


V'V  +  V-W-v1  =  0 


'  I  lio) :  I  '/*  1  (<0) 


const . 


(3-16) 


2060 
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W  +  l7K(a0,illl 


2(:0>’  1  '7  1  I  (;0>  l’l  Ct 


B,(a0)  +  l',  fA]  (a0)a|)  +  *-’7  fA2<a0)a  1  1 


(3-1 7) 


f 3-18) 


^1‘V  - ijfvv  =  r!  =  c°nst-  >■  •  018) 

l(V  =  cl  =  conSt‘ 

System  (3-16)  or  system  (3-17)  ran  be  rejected  if  its  solution  is  such  that  the  expressions 
'  <  ‘  )  :  and  )f  have  equal  orders. 


Let  us  suppose  that  all  values  of  ,  $  are  found  up  to  the  number  n  =  k  -  I 


We  write 


k-  i 

y  -  > 

Li  n  n 


(3-19) 


Then,  z  =  /  +  z  . 


Repeating  the  calculation,  we  obtain: 


+  VWV  -  0 


...  .  /  ■  .  k-  I  - 

'  I  '  '  1 1 1  '  k  :/  'k  ' 


=  const. 


(3-20) 


>  +  '  T  !A',(an)a  j  -  0 


.  ,  ,  k-i. 

•  ,<  i/  ..  <  ‘  )  =  c.  =  const. 

0  k  /  k  k 


vjk_l)  +  vwv  +  t-k  -Aj(ao)ak^  =  0 
^vv/y;H>  =  ck  =  const 


(3-21) 


(3-22) 


l-:2(*0)V/;k<t  >  =  Ck  =  COnSt’  J 

k-  I  k-  I 

where  a  and  f  are  vectors  with  the  coordinates  a  ,  -fn  ;  0  <  n  <  k  -  I  ; 

B,  (a^  7.  (  ; /-  )  is  the  dominant  member  of  the  expression  F  z7  *  . 

k  k 

Systems  (3-lb)  and  (3-17)  can  be  rejected  on  the  basis  of  the  same  considerations. 


..  '■  :  The  cxpans  i  on  (1-2)  under  E  he  font]  i  l  i  on  s  of  (  1-4)  and  lor  /.iron  .nines  ui 

is  unique  ,  ’■  i'  two  functions  coincide  to  the  desired  extent  for  tin-  sir, ill  t mite 

n 

0  •>  s  >0  ,  then  both  they  and  their  exp.msions  coincide  over  the  whole  range  ot  the 

de t  erminnt ion . 

It  is  sufficient,  therefore,  that  the  properties  of  sub-points  (I)  to  II)  should  he 
maintained  merely  tor  small  values  of  s  . 

The  properties  of  sub-points  (2)  and  (1)  are  directly  associated  with  the 
existence  of  the  Frdchel  derivative  of  the  operator  F  .  If,  for  example: 

Fix'  +  z")  =  Fz'  e  r'(x')z"  +  o(z")  ,  (  ,-j'l) 


where  F’(z')  is  ;t  linear  uniform  operator  depending  on  z  (Freehol  derivative;,  then 
the  properties  of  (1-7)  and  (2-8)  are  maintained  if: 

1  innF'  (z  |  )z"  \  j  {?'  (?..',  )z"j  =  I  for  z|/zj  -■  I  ,  s  •  0  .  <’  1-24) 

Thereupon  the  following  formulae  are  valid: 


DP  Fz  =  Df>  F(a  $0)  , 


0-25) 


DP  (Fz  -  F(z  -  7^)1 


DP  f  F 1 


(ao;o)ak  'V 


0-26) 


4  Applying  the  above  considerations  to  equation  (2-10),  let  us  first  of  all  formulate 
additional  restrictions  of  a  physical  character: 

(i)  the  absence  of  self-intersection  of  the  tangential  discontinuity: 


(2) 


z  (-. 

1  ,s)  -  Z (<i'0 , s) 

i/l*(i,s)|  >  f(., 

14-1) 

noU’: 

rn  no  r 

t* 

tionr  hv'f  i  ■  •  >  ■ 

•.  ■  (■-:)  v 

J  (■-!).) 

when*  )<•:.. 

r'V 

is  a  certain 

function ; 

f  F, 

1  ,v2) 

at  0 

for 

1  *  -2  ; 

the  boundedness  of  the  velocity  jump  at  the  tangential  discontinuity: 


a  z  ( :,. ,  s )  /  ti. 
z ( 1 , s) 


>  c  '  •  0 


where  c '  is  a  certain  constant; 

(5)  the  boundedness  of  the  total  circulation: 


(4-2) 


0-3) 


(4)  the  tangential  discontinuity  springs  from  t  lie  edge: 
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1  I 


2  ( 0 ,  s )  =  0  ; 


( 4-4) 


t5)  moreover,  we  suppose  chat  the  series: 


Ed  a 

n=0 


da 

n  )7- 

df,;  *n  3u 


(4-5) 


converges  uniformly  with  respect  to  u  . 

From  (4-1),  (4-2)  and  (4-5)  it  follows  that  for  any  values  of  .  , 


Z|  (<jj(  s)  -  z  (<^2s) 
lim  - t — r - - - — v-  =  0 

Vo(“i)  "  ♦oW 


for  s  -  0  .  (4-6) 


Therefore,  for  small  values  of  s  : 


2(00., s)  -  z(u>,s) 


“  a0(u)^ 


»  z  (w*, s)  -  z } (w, s) 

\  «0^a0(“*)  "  a0(u)))  / 


*o^ao(“*)  '  ao(li,)) 


1  + 


/  z| (u*»s)  -  Z) (u,s)\k 


(4-7) 


This  makes  it  possible  to  write  p(z)  in  the  form  (3-6).  An  analogous  result  is 
obtained  for  the  operator  E  ,  and,  consequently,  for  all  the  left-hand  side  of  equation 
(2-10).  The  properties  (3-7),  (3-8)  of  the  left-hand  side  of  (2-10),  written  in  the  form 
(3-6)  can  be  directly  verified.  Here  formulae  (3-25)  and  (3-26)  may  be  used.  Denoting 
by  the  symbol  F  the  operator  included  in  the  left-hand  side  of  (2-10),  we  obtain: 


As  in  the  case  of  point  (3),  we  obtain  three  equations,  two  of  which  can  be  rejected  not 
only  on  grounds  of  the  conditions  of  uniqueness.  For  example,  let  us  look  at  the 
equation: 


(*' 


0 


(4-9) 


By  virtue  of  (4-3)  and  (2-11), 


V 


12 


j 

-I 


* 


di.i 


a0* 


*  0  , 


therefore  ( 1  -  u)  a.  +  a„  =  0  ,  a  =  c , ( I  -  . )  ,  c ,  =  const.  ; 

du>  0  0  o  !  I 

which  contradicts  condition  (4-4).  (FuH  tor  'o  not-  :  w  n»'  nor  -v<  rr-ationr  here.) 


The  second  equation: 


(4-10) 


has  the  following  solution: 


which  does  not  satisfy  condition  (4-2),  since  aa^/Sw  =  0  for  w  =  I  .  In  this  way,  if 
a  solution  of  the  problem  in  the  form  (3-2)  exists,  then  satisfies  the  system 

(3-12).  In  so  far  as  4> j 4>q  =  4)^  (editor's  note:  the  original  has  <J>a  J,  and 

<t>2<t>0  =  ’  then’  for  *0  we  obtain: 


(4-11) 


where  4>q  =  di^/ds  . 

Since  the  constant  factor  associated  with  4>n  is  insignificant,  then,  for 

0  ^ 

definiteness,  it  is  possible  to  take  c ^  =  j  ,  which  by  virtue  of  (3-12)  unambiguously 
determines  the  equation  for  a^  . 

Transferring  attention  to  the  desired  ^  .  . . ,  4^  we  first  obtain: 


Dp(fz  -  F(z  -  zk)} 


(4-12) 
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t  ,  then  with  an  acruraev  of  up  to  the  value  ol  a  constant  multiplier  of  the 
S  t  nee  :  =  s  ' 

2  .4 

function  -q'q  *s  equal  to  the  function  1/ : ^  .  Theretore  we  have  the  situation 
described  in  point  (3).  Thus: 


1  (if  )  f 
I  yv0  k 


Extracting  the  dominant  member  from  the  expression  Ff  1^:^)  >  bearing  in  mind  equation 
(3-12)  and  the  function  =  s 3  ,  we  obtain: 


i 

DP  F(aQs  T] 


_  t_ 

Thus,  y j  v  q  =  s  ’  .  Investigating  the  systems  (3-16)  to  (3-18)  according  to 
principles  laid  down  in  point  (3),  the  first  two  systems  can  be  rejected.  From  the  third 
system,  with  an  accuracy  of  up  to  the  value  of  a  constant  multiplier,  we  obtain: 


f> |  ( s)  =  s  .  (4-13) 

Applying  the  method  of  mathematical  induction  for  ,  we  find: 

bk(s)  =  s(2k+l)/3  k  =  0,1,2 .  (4-14) 

Introducing  the  result  found  into  (2-11): 


no 


I 

n=0 


(2n+l )/3 

y  s 


Y  =  const, 
n 


(4-15) 


From  (2-12),  it  follows  that: 


oo 

'.E 


2n/3r 


AD  =  pv  )  s  D  ,  D  ■  const. 
"  n  n 

n=  1 


(4-16) 


Using  (4-16),  (1-2),  (1-5),  (1-6)  and  (2-7),  we  obtain 

CO 

Ev 

n*  1 

00 

a>  Y  ■  ■  D  (2n 

/  ,  2n  +  3  n 


c  =  aAu/2  +  JaA 

y 


(2a/A) 


2n/3 


/A) 


2n/3 


n*  l 


(4-17) 
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For  non-steady  flows,  we  use  (1-3)  and  (1-4);  accordingly,  we  obtain  c  ,  >1  and 
M"  . 


0 


+  ']Dn(2 


2a  F../>) 


2  n  /  3 


M*  * 

7 


0  D  (2 af/>.) 
n 


2n/3 


(4-18) 


M" 

z 


a> ~/2  +  Ja.l  ^  Dn(2ar/t)2n/3 
n=  I 

aAf'  v2aC/A)2n/3  ,  t  =  t0v.,/h0 

n=  I 


(4-19) 


In  order  to  find  the  coefficients  ,  it  is  necessary  to  solve  the  equations  for 

a  in  systems  (3-12)  and  (3-22).  However,  since,  when  one  uses  numerical  methods  with  the 
n  5 

application  of  the  regulartsat ion  principle  ,  there  occurs  a  convergence  not  only  in  the 

case  of  the  desired  functions,  but  also  in  the  case  of  all  their  derivatives,  then  the 

coefficients  D  can  be  found  by  means  of  a  treatment  of  the  results  of  the  numerical 
n 

calculation  of  the  initial  problem  ((2-1)  to  (2-5)).  In  Fig  I  is  shown  the  dependence, 
found  by  the  numerical  calculation  of  the  initial  problem,  of  Ac  rx  ■’  A  3  as  a  function 
of  2aA  . 

It  is  noticeable  that,  at  point  0  ,  there  is  a  vertical  tangent.  This  same 
quantity,  as  a  function  of  (2a/A)^  does  not  have  a  vertical  tangent  and,  in  the  interval 
(0,1),  it  can,  with  a  high  degree  of  accuracy,  be  found  approximately  using  the  following 
polynomial  (cf  Fig  2): 

- *-  =  3.61  +  1.206  -  0.556  -  0.256  ,  0  =  .  (4-20) 

it  A 

a  A 

The  coefficients  of  this  polynomial  determine  the  first  four  coeficients  : 

D|  =  4.55  ,  =  1.56  ,  =  -0.692  ,  =  0.314  .  The  constant  c^  in  the  formula  of 

A. A.  Nikol’skii  is  numerically  equal  to  the  first  coefficient  of  the  polynomial  (4-20). 

A  comparison  of  results  of  the  calculation  and  of  experiment  is  given  in  Ref  5. 

In  conclusion,  two  observations  may  be  made. 

Equation  (4-10)  corresponds  to  the  case  when,  in  the  flow,  the  pressure  is  admitted 
of  a  point  force  applied  to  the  end  of  the  tangential  discontinuity.  Thereupon,  the 
function  can  be  determined  if  one  defines  the  law  of  change  of  this  force. 
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Hy  virtue  of  formulae  (4-1/),  there  is  the  following  diiterenli.il  re  I  .U  i  nnsh  i  p 

hi' tween  M  and  e  : 
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